arXiv:1502.02820vl [hep-th] 10 Feb 2015 


KEK-TH-1794 

AP-GR-119 

OCU-PHYS-419 


Quasinormal modes of (Anti-)de Sitter black holes 

in the XjD expansion 

Roberto Emparan“’^, Ryotaku Suzuki*^, Kentaro Tanabe'^ 

‘^Institucio Catalana de Recerca i Estudis Avangats (ICREA) 

Passeig Lluis Companys 23, E-08010 Barcelona, Spain 
^Departament de Eisica Eonamental, Institut de Ciencies del Cosmos, 
Universitat de Barcelona, Marti 1 Eranques 1, E-08028 Barcelona, Spain 
'^Department of Physics, Osaka City University, Osaka 558-8585, Japan 
‘^Theory Center, Institute of Particles and Nuclear Studies, KEK, 

Tsukuba, Ibaraki, 305-0801, Japan 

emparcLn@ub.edu, ryotaku@sci.osaka-cu.ac.jp, ktcLnabe@post.kek.jp 


Abstract 

We use the inverse-dimensional expansion to compute analytically the frequencies of a set 
of quasinormal modes of static black holes of Einstein-(Anti-)de Sitter gravity, including 
the cases of spherical, planar or hyperbolic horizons. The modes we study are decoupled 
modes localized in the near-horizon region, which are the ones that capture physics peculiar 
to each black hole (such as their instabilities), and which in large black holes contain 
hydrodynamic behavior. Our results also give the unstable Gregory-Laflamme frequencies 
of Ricci-flat black branes to two orders higher in 1/D than previous calculations. We 
discuss the limits on the accuracy of these results due to the asymptotic but not convergent 
character of the 1/D expansion, which is due to the violation of the decoupling condition 
at finite D. Einally, we compare the frequencies for AdS black branes to calculations in 
the hydrodynamic expansion in powers of the momentum k. Our results extend up to 
for the sound mode and to k^ for the shear mode. 



1 Introduction 


The quasinormal modes of a black hole spacetime encode important aspects of its dynam¬ 
ics [1]. In particular, the AdS/CFT correspondence implies that the quasinormal modes 
of Anti-de Sitter black holes describe the relaxation to thermal equilibrium in the dual 
field theory [2], As argued in [3] and further elaborated in this article, the limit of large 
number of dimensions D isolates a subset of quasinormal modes associated to particu¬ 
larly interesting black hole dynamics, and allows efficient analytic computation of their 
frequencies. 

In more detail, when D is very large the gravitational field of a black hole gets strongly 
localized close to the horizon [4, 5]. The existence of a well-defined ‘near-horizon region’ 
[6] splits the quasinormal spectrum into two distinct sets that capture very different black 
hole dynamics. The most numerous set (oc D^) are ‘non-decoupling’ modes that straddle 
between the near-horizon zone and the asymptotic zone; they are largely insensitive to 
the peculiarities of the black hole, which for these modes is simply a hole in a background 
spacetime [7]. In contrast, the much smaller set (oc D^) of ‘decoupled’ quasinormal modes 
have support only in the near-horizon region, where they are normalizable states to all 
orders in l/D] they capture properties specific to each black hole, such as the instabilities 
of certain higher-dimensional black holes and black branes, and the hydrodynamic modes 
of black branes [5, 8].^ 

Due to the simple form of the near-horizon geometry, the decoupled quasinormal fre¬ 
quencies can be calculated perturbatively in analytic form to several orders in 1/D. Fur¬ 
thermore, the universality of the leading-order near-horizon geometry implies that the 
structure of the calculation is essentially the same for all static neutral black holes, be 
they asymptotic to de Sitter, Minkowski, or Anti-de Sitter, with spherical, flat or hyper¬ 
bolic horizons. This allows us to obtain unified analytical formulas for their frequencies, 
valid up to next-to-next-to-next-to-leading order (3NLO) in 1/D, and for planar horizons 
up to 4NLO. Via the AdS/Ricci-flat correspondence [10], the latter also give the unstable 
Gregory-Laflamme (GL) frequencies of Ricci-flat black branes [11] to two higher orders in 
1/D than in [5]. 

For AdS black branes a hydrodynamic gradient expansion has been applied to obtain 
the sound- and shear-mode dispersion relations uj{k) analytically (exactly in D) up to 
for the sound mode, and up to k^ for the shear mode [12]. Where there is overlap, our 
calculations perfectly match these results. As we will see, the regime of applicability for 
the large D expansion is not essentially different than for the hydrodynamic expansion. 
However, the large D expansion allows a much simpler calculation of higher orders in k: 
our results extend up to k^ for the sound mode, and up to fc® for the shear mode. 

It is natural to ask whether these higher-order corrections always increase the accuracy 
of the results when finite, possibly low, values of D are substituted. In other words, is 

'^The decoupled spectrum at large D was first identified numerically in [9]. 
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the l/D expansion for the decoupled quasinormal spectrum a convergent series? We will 
explain that it is not, bnt is instead only asymptotic. The decoupling condition, which 
holds to all pertnrbative orders, is violated by effects that are not perturbative in l/D. 
Nevertheless, we obtain excellent agreement for certain quantities, snch as the GL critical 
wavelength, even down to D = 6, which we reproduce within ~ 2% accuracy. 

In the next section we introduce the class of black holes we study and their large D 
limit. In sec. 3 we explain qualitatively the appearance of the decoupled spectrum and 
its main features. The main resnlts, namely the scalar and vector decoupled frequency 
spectra and the GL unstable frequencies, are given in sec. 4. Sec. 5 examines the lack 
of convergence of the 1/D series for deconpled modes. Sec. 6 compares onr results with 
those obtained exactly in D in the hydrodynamic expansion at low momenta. We conclude 
in sec. 7. Appendices A to D contain further technical details and lengthy results. In 
particular, appendix D gives explicit expressions for the frequencies for spherical horizons. 


2 Set up 

We consider (A)dS black holes with metric 

ds"^ = -f{r)dt^ + f{r)~^dr^ + dcr|-,D_ 2 , 

with 


/(r) = iL-Al2-(^ 


D-3 


( 2 . 1 ) 
( 2 . 2 ) 

{D-l){D-2) 

parametrizes the cosmological constant A in units of the black hole size-scale tq. The 
solution has another, discrete parameter. 


and where 


A = 


2Arg 


K = ±1,0 


(2.4) 


for the cnrvature of the metric da\ j ^_2 on the ‘orbital’ space ICd- 2 '- a unit sphere, a 
plane, or a hyperboloid, for K = ±1,0, —1 respectively. 

It is convenient to use 

n = D-3 (2.5) 


instead of D as the perturbation parameter, and 



( 2 . 6 ) 


as the radial variable appropriate for the near-horizon region, which is defined by R ^ e"". 
In terms of it, the metric (2.1) is 


ds^ 


-/(R)dt2 + !±r2A 


dR2 

wm 


± R2/”rg da 


2 

K,D-2 


(2.7) 


2 




with 


/(R) = K- 


( 2 . 8 ) 


/(R) = K- AR2/- - 1. 

The horizon at R = Rh such that /(Rh) = 0, has surface gravity k given by 


nK 


nro = 


2R 


1/n 

H 


-(f + l)ARr 


(2.9) 


Taking n —)• oo with R and A finite, to leading order the geometry in the (t, R) directions 


is^ 


with 


ds^2d) = -i —-7:^]dt + 


Rn R 


dR2 




R2 


1 


Rn = 


1 


( 2 . 10 ) 


( 2 . 11 ) 


K-X 

If we rescale {t, R) —)• (Rot, RqR) the metric becomes the same for all values of K and A. 
However, subleading corrections distinguish among them, as we can see in the horizon 
location which is at 

Rh = Rq T + O (n ) . (2.12) 


n 

In the planar case = 0 all values of A < 0 are equivalent to A = — 1. In the hyperbolic 
case {K = —1, A < 0), when n is finite there exist black hole solutions with negative mass 
parameter rg < 0, with the lowest negative mass corresponding to an extremal black hole. 
It is easy to see, however, that as n —)• oo the range of allowed negative masses shrinks to 
zero, and as a consequence A is bounded above by A < —1. More generally, the large n 
limit yields a good near-horizon geometry only when 


Ro>0. (2.13) 

Note that the hyperbolic slicing of AdS (obtained for rg = 0) lies outside our study, which 
is expected since there is no strong localization of the field close to an acceleration horizon. 


3 Decoupled spectrum: qualitative aspects 

Before entering the calculational details, let us explain the existence and basic properties 
of the decoupled quasinormal spectrum at large D. 

At first sight, one would not expect any decoupled dynamics in the near-horizon region 
of neutral static black holes: this region has only a very short radial extent ~ vq/D away 
from the horizon, in contrast to the long throats of (near-)extremal black holes which 
support excitations that can propagate inside them without ever leaving. Indeed, the 
excitations with frequencies ui^h ^ D/tq characteristic of the near zone are not decoupled: 
they are non-normalizable states of the near-horizon geometry [3]. However, as we will see 

^We omit here the ICd- 2 , whose size plays the role of the dilaton for the effective two-dimensional 
gravity [6]. 
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Figure 1: Radial potentials F(r») of gravitational vector (left) and scalar (right) perturbations, 
for n = 500 and harmonic eigenvalue = 2n, for representative examples of (A)dS black holes. 
The plots show the potential around the region where the near and far zones meet (the horizon is 
at r* —^ —oo). Decoupled quasinormal modes are localized around the minima in the near-zone. 
The potentials for tensor perturbations differ only in the absence of these minima. 

presently, this geometry supports a few modes that, to leading order, are not dynamical'. 
in near-horizon scales they are zero-modes with frequency ~ —)• 0. Furthermore, 

the universality of the leading-order near-horizon geometry implies that these modes are 
present in all neutral black holes. 

Such leading-order zero-modes exist for gravitational perturbations which are vectors 
or scalars of the space JCd- 2 , and are absent for gravitational tensors (and also for free 
massless scalar fields, which satisfy the same equation). To see this point, note that the 
wave equation can in all cases be written in the form [13] 

= 0 (3.1) 

where = f dr/f{r) is the conventional tortoise coordinate and 

s = 0,l,2 (3.2) 

denote gravitational scalars, vectors and tensors, respectively. The explicit form of the 
potentials Vg will be given later below; their form at large D is presented in fig. 1 for 
vectors and scalars, zooming in around the most interesting region, r* = 0(1), where 
the near and far zones meet. At large r* the potentials reflect the properties of the 
background: decaying for asymptotically dS and flat space, and growing (‘bounding box’) 
for asymptotically AdS. In contrast, near the horizon the form of the potential becomes the 
same independently of the cosmological constant or the horizon curvature. In the direction 
of decreasing r*, as we enter the near-horizon zone the potential drops exponentially from 
a height ~ /rq. In this near-horizon region the potentials for gravitational scalars and 
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vectors have negative minima, which are absent for tensors.^ These minima allow the 
existence of non-trivial solutions to the equation (3.1) with a; = 0. Since at large D both 
d?/dr1 and 14 scale like D'^jrQ, these zero-modes are actually leading-order solutions to 
(3.1) for modes with frequency u = o{D/ tq). 

These leading-order static modes become dynamical, i.e., their frequency is non-zero 
uj ~ 1/^0, at the next order in the 1/D expansion. At this order, the differences between 
black holes (due to the cosmological constant or the horizon curvature) make themselves 
present and therefore the frequencies differ for each of them. However, the structure of the 
perturbation expansion is the same for all of them since it is controlled by the leading-order 
solution. 


4 Decoupled spectrum: analytic results 

Henceforth in this section we use units where 


ro = 1. 

The explicit form of the potentials 14 is [13] 

/(R) 


v :,=2 = 


R2/r, 


q^ + 2K + 


R- 


dR 


-b 


14=1 = 


/(R) / 2 3 -b 1)^ 1 — 1 

' ' -K — 


r2/. 


and 


q^ + 


14=0 = 


R 



(4.1) 

n? — 1 , \ 


4 

(4.2) 


(4.3) 


/(R) 


Q(R) 


R2/. 




2 ’ 


where m = q^ — {D — 2)K and 
(n -b l)^(n -b 2)^ 


Q(R) = 


-b 


R3 

(n -b l)(n -b 2) 
R2 


(4m(2n^ -b n -b 3) -b (n -b 2)(n — 3)(n^ — 1)K) 


12(n + l)m ^ ^ 


R 

-b4(n -b 3)(n -b l)Km^ + 16m' 

\2t^ , 113 


(4.4) 


(4.5) 


(<":±.?1(?±|),(»±0. _ + 4(„ - 1)(,. - 3)md ARV» . 


We use q^ to denote the eigenvalues of the Laplacian on K,n+i- It is easy to see that 
the minima of the potentials and hence the decoupled modes only exist for s = 0,1 and 

®Fig. 1 leads one to expect that the 1/n expansion is more accurate for vectors than for scalars, which 
indeed appears to be the case. 
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for eigenvalues such that Jit? 
renormalized eigenvalues 


In particular, 


0 [3]. Thus we take = 0{n) and introduce finite 



(4.6) 


+ (^ = 1 ) 

f (^ = 0) 


(4.7) 


where i = 0{vP) is the angular momentum number on and k = 0{^/n) is the 

momentum along For K = —1, if the hyperboloid is not compactified the spectrum 

is continuous like in the planar case, but we will not be more specific about this. 

Quasinormal modes are solutions to (3.1) that satisfy specific boundary conditions. 
For decoupled modes, these are as follows. At large distance in the near-horizon geometry, 
where 1 ^ R <C e", we impose 


Ts(R —y cx)) —)■ 


1 


(4.8) 


at all orders in the 1/n expansion. This expresses that the mode is normalizable in this 
geometry (non-normalizable modes are oc \/R at large R). Since it is (perturbatively) 
shielded from the far zone, it can be matched to a purely outgoing wave in a Minkowski 
background [3] or to other suitable function for other asymptotics. 

The solution must be ingoing at the future horizon. This is achieved when 

T,(R) = (R - R^)-*-/(2-) ^^(r) (4,9) 


where (^s(R) is regular at R = Rh- 
Observe that 


• the condition at R —?• oo is independent of the asymptotics in the far zone; 


• since k (x n, the horizon condition to leading order in 1/n is simply finiteness of 
^s(Ro). 


These two properties, together with the universal character of the geometry (2.10), imply 
that the leading order solutions (the zero modes) are the same for all values of K and A, 
which allows to study these quasinormal modes in a unified manner for all the metrics 
( 2 . 1 ). 

The perturbative calculation of the decoupled quasinormal frequencies now proceeds 
as for Schwarzschild black holes {K = 1, A = 0, Rq = 1) in ref. [3]. Like in that instance, 
for the scalar perturbations it is more convenient to use a formulation in other variables 
than 4 's=o(R)) which we give in appendix A. We omit the lengthy but straightforward 
details of the calculations and just give the final results. Writing 


U} = 


E 


n* 


(4.10) 


we obtain 
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Vector-type quasinormal frequencies: 

^{0) = i{K - q^) , (4.11) 

W(i) = -i{K - f) (In Ro + 2), (4.12) 

^(2) ^ ~ (61nRo + 7r^) - - 31nRo (In Rq + 4) + 2 (tt^ - 12)) ,(4.13) 

<^(3) ^ ~ Ro + 6) - 27r^ -h 24) In Rq - 8 (tt^ - 6) 

+2Ro (^12RoC(3) (A - f) {K - f) + (In Rq + 4 - 2ARo (In Rq + 2)) 

-|-2A^Ro (In Rq + 2) (3 In Rq + vr^) + (vr^ — 12) A In Rq^ • (4.14) 


Vector modes are purely imaginary. For K = 1 a. perturbation with (f = 1, i.e., f = 1, is 
not a quasinormal mode but corresponds to adding angular momentum to the black hole. 

Scalar-type quasinormal frequencies: 




Ret(;(i) = 


{K-e) (iF(lnRo + 2)-a-^ 


X-iF2 


(4.15) 


(4.16) 


Rew(2) = -- 




^ 8R:3Ro (^k - f) {2K (3 In Rq + vr^) - 7r^f) 

-4 ^3R:^ (In Ro (In Ro + 8) + 8) + 6K^f (-2 In Ro + vr^ - 7) 
—K'^q^ (3 In Rq (In Ro -|- 8) -|- 87r^ -|- 3) 

+2Kq^ (3 In Rq -|- vr^ -|- 9) — 3g®^ 

+ ^(3K^ (In Ro (2 In Ro + 11) + 11) 

Ro V 


(4.17) 


+Kf (-3 In Ro (2 In Ro 17) -h 27r2 - 69) 


+q^ (l2 In Ro — 27r^ -|- 


Rg ’ 


lmuj^Q)=K-f, (4.18) 

Imwq) = (In Ro -|- 3) — RT (In Ro -|- 2) , (4.19) 

Im W( 2 ) = -y (iF - f) - 2Rr (3 In Ro vr^)) -k y (In Rq (In Ro + 8) -k 8) 


-y (In Ro (In Ro -k 10) -k 14) . 


(4.20) 
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The expressions for cn( 3 ) are very lengthy and we defer them to appendix C. The K = I 
scalar modes with i = 1 are gauge modes. 

These results are valid for all K and A satisfying (2.13). For ease of reference, we give 
the expressions for the particnlar case of K = 1, i.e., (A)dS Schwarzschild black holes in 
appendix D. For A = 0 Schwarzschild black holes, these results reproduce those of [3]. 

Large black hole limit vs. large n limit. At any finite n, the large black hole limit 
—A = —)• oo of the AdS Schwarzschild solution, with K = 1, resnlts in the AdS 

black brane. However, when we compute modes of AdS Schwarzschild black holes in the 
1/n expansion we assume that tq/L remains finite as n —>• oo (more precisely, we require 
that rg/L^ <C e”). Thus the two limits n —?• oo and ro/T —?■ oo need not commnte. Indeed, 
while the leading order terms a;(o) for AdS black branes are correctly obtained as the limit 
vq/L —>■ oo of AdS Schwarzschild frequencies, the next-to-leading order corrections are 
not. 


4.1 AdS black branes 

For planar horizons it is possible to perform the analysis in a different manner [14] (see 
appendix B): given a momentum vector k°‘ along one can decompose the perturba¬ 

tions into scalars, vectors and tensors of its little group SO{n). This has allowed us to 
carry out the expansion to one higher order than the previous results, up to 1/n^. 

For this case, the expansion parameter is more appropriately 1/(11 — 1) = l/(n -|- 2) 
rather than 1/n. So here we introduce 


n = D — l = n + 2, 


(4.21) 


and 


P = ^ = ( 1 - T 


n 


n 


and express results as an expansion in 1/n. We find: 


(4.22) 


Vector (shear) mode: 


Lo = —ik^ 


1 + 


PP 4P(1 + P)C(3) , 4PP(1 + 7P + P) , 

T j — A T d(n ) 


3n2 




45P 


1-2 f.2 I f.4 12 I 7p I f.6 

= -iP 1 + 2C(2)^ - 4C(3)P^ + 8C(4) - .— _ + ^ ^ + ©(h'^) 


n^ 


n'' 


n’ 


(4.23) 











Scalar (sound) mode: 


Rew± = ±/c 


1 + 2P 1 


2n 


1 


\ 16 


1 / 35 


+ ^2 0-^ + 

\ 8 


9 vr 




+ 3^ I ~ — k I “ H—~—h 4^(3) ] k I T + ~ 2(^(3) j H—— 


8 6 


25 37r^ 47r^ 


k^ 


+ + ^ (-T1J- —+ ^ + 2C(3)) +A: (-- —+ 


16 


45 


13 3tt^ 


16 


297r'* 


-5C(3) 


57r^ vr' 


+ A) I —— --—h — — 22(^(3) j -— 1 + 0{n 


15 


5^8 




(4.24) 


and 

Imcjj- = —P 
+ 


1--+4 

n 


2^2 


-1 + 


TT^k 


+ 


n"’ 


-1-k^ 




+ 8C(3) - 4fc4C(3) 




-1 -r 


- 16C(3) ) + 


3l7r^ 


+ 36C(3) + 


47r^fc®\ 


45 


- 


+ 0(n-5 

(4.25) 


The appearance of the ^ function in these series will be explained in sec. 6. 


4.2 Gregory-Laflamme unstable frequencies 

The AdS/Ricci-flat correspondence of [10] relates the quasinormal spectrum of H = h + 1 
AdS black branes to the spectrum of fluctuations of Ricci-flat black p-branes in dimension 


D = n + p + 3 . 


(4.26) 


The scalar sector of the latter is known to contain unstable modes [11]. 

According to [10], the map requires replacing h —?• —h, so we also have k —>■ ik. By 
applying this to eqs. (4.24), (4.25), we find the imaginary frequencies as 

k / ... 

= zbfc — k"^ -(il -|- 2/c =p 2k^ 

^ 2n\ 


+ —r^ (±9 + 24k ± 12P =F Stt^P Sir^P =F 12fc‘^ 
24h2 V 




k 

4h 

k 

48n^ - 


=F15 - 48fe =F 2(27 -b dvr^ -b 96C(3))P -b 64(7r2 -b 6C(3))P 
T12(3 -b 4p - 8C(3))P - 192C(3)P ± 24pT 


~*~n^ 


TT^ vr^ 


+ + 2C{3)) +k^-- + - + 16C(3) 

128 V 16 8 45 ^ \ 3 9 ^ 


-4 , 13 37r2 29P A /3lP 


, 5 57r^ P \ 4pP 5P 

±P - -b--b 22C(3) H-=F- 

4 6 15 '’V 7 ) w -r g 


+ 0{n 




(4.27) 
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0+ is the unstable mode on the black brane. We have also obtained this same result by 
solving directly the perturbations of the Ricci-flat black brane, thus managing to extend 
the calculation in [5] to two higher orders. 

Using (4.27) we can find the critical wavenumber at the threshold of the instability, 
i.e., /cgl = such that U+(/cgl) = 0{n~^). We find 

*=0. = vs (l - i ^ + (2C(3) - I) ^ - 5C(3)) ^ + 0(«-^)) . (4.28) 

In fig. 2 we compare this result to the values found in [4] from the numerical solution of 
the problem. For n = 2, e.g., a black string in D = 6, the numerical value is kcL = 1.269, 
while (4.28) gives = 1.238, which is off by 2.4%. 


I<GL 



Figure 2: Wavenumber kcL of the marginally stable mode of black p-branes as a function of 
n = D — p — 3. We plot successive approximations to eq. (4.28): to 1/n (purple dotted line); 
to 1/n^ (red dash-dotted line); to 1/ri^ (blue dashed line); to 1/# (black solid line). Dots are 
numerical values from [4]. Units are rg = 1. 


5 Non-perturbative breakdown of decoupling 

We now argue that the l/D expansion for decoupling modes does not converge but is 
instead only asymptotic, since effects non-perturbative in 1/77 arise from the breakdown 
of decoupling at finite D. 

Considering for simplicity static perturbations of a Ricci-flat black brane, the solution 
in the far zone is, up to normalization, given by a modified Bessel function, 

^'(r) = \/rKfi +2 (kr). (5.1) 

2 

The Bessel function at large order has the asymptotic expansion 

Ky{x) ^ (i/2> 1) (5.2) 

i>0 
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where Oj is a polynomial in x, behaving like a* ~ x* so in the overlap zone we have, 
schematically, 


'l'(r) ~ 



(5.3) 


We see that 


Zero <C n , 


(5.4) 


is required for 'k(r) to satisfy the condition (4.8) that the mode decouples. This is a limit 
on the range of momenta to which the large h expansion for black branes is applicable, 
and it is valid at all orders in the 1 /h expansion. This bound is indeed consistent with the 
fact that in eqs. (4.23), (4.24), (4.25) each additional order in 1/h brings in an additional 
power e. 

However, large orders in the expansion (5.3), with i = h+0(l), give a behavior T rsj Vr 
which violates the decoupling condition (4.8). This breakdown is non-perturbative in 1/h, 
and gives a value of 


non-perturbative corrections 



(5.5) 


Generically, for planar horizons (Ricci-flat and AdS black branes) we have Zero ~ v/h so 
we expect 

non-perturbative corrections = O {K = 0), (5-6) 

This implies that, for h sufficiently large, the perturbative expansion is reliable only for 
orders m such that 

71 

m<- + 0{l). (5.7) 

This rough estimate agrees fairly well with the fact (not shown here) that the inclusion 
of the 1/h^ correction in eq. (4.27) gives a better overall fit to numerical calculations of 
the curve Q+{k) only when h > 8.^ Note however that (5.5) implies that the accuracy at 
small k can be considerably better than indicated by (5.7). We elaborate this point in the 
next section. 

For spherical black holes, a similar argument applies replacing Zero —S’ ^ = 0(1), so the 
expected corrections in this case are 


non-perturbative corrections = O (n ”) (K = 1). (5.8) 

These are smaller than (5.6), which agrees well with the good accuracy found at low n in 
[3] for the computation of decoupled quasinormal frequencies of Schwarzschild black holes. 

To finish, let us mention that the numerical calculations for AdS black branes in [14] 
of the sound (scalar) frequency in Zl = 5 suggest that at large momenta. Re —>■ k and 

'^However, the best approximation to koL (although not to fi+(k) overall) for h = 2,3 is the fourth-order 
one, as fig. 2 shows. This might be a fortunate accident at low n, since for 4 < h < 8 the second-order 
result for koL is (slightly) better than the fourth-order one. The third-order correction does not improve 
on the second-order result until n > 12. 
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that |Imw±| decreases. This is not visible in our results, so presumably these effects are 
non-perturbative in 1/h. 


6 Comparison to hydrodynamics 


Although the 1/D series of quasinormal frequencies is not convergent, certain terms in 
them can be resummed to all orders in 1/D. This is the case at least for vector and scalar 
frequencies when expanded at low momenta. Indeed these terms have been computed 
exactly in D in the hydrodynamic expansion [12].^ The vector (shear) mode frequency for 
an AdS black brane in D = h + 1 dimensions (h > 3) is® 

p / jU2 \ 

LO = —i— ( 1 + V2-f^2/fi ) + 0{k^) (6.1) 

n \ ' J 

where Hm is the m-th ‘harmonic number’. The scalar (sound) mode frequency is 


U!± = ±- 


k 


n — 2 


— i- 


,2 + -^a/n) , 3 


■P ± 


yjn — 1 n{n — 1) n^{n — l)^D 

These results can be expanded at large n using 


r + 0(P). 


( 6 . 2 ) 


D2/n — 



(6.3) 


which is a convergent series for |h| > 2, and then we find perfect agreement with our large 
n calculations in eqs. (4.23), (4.24), (4.25). It is amusing to note that the appearance of 
the function i^ 2 /n these dispersion relations could have been obtained by resumming 
(using (6.3)) the natural guess for the entire l/fi series from the terms present in (4.23).^ 
Since the temperature of the black brane (with rg = 1) is T oc h, the range of applica¬ 
bility of the hydrodynamic expansion, namely fe <C T, coincides parametrically with the 
range of momenta to which the large h results apply (5.4). In this view, the large h ex¬ 
pansion does not afford a larger regime of applicability than the hydrodynamic expansion. 
The main advantage of the large h expansion is that it allows a simpler computation of 
higher powers of the momenta (in our case, up to k^ for the sound mode, and k^ for the 
shear mode), some of which can be relevant (barring non-perturbative effects) at relatively 
low dimensions. These results also contain information about the transport coefficients at 
higher orders in the gradient expansion. 


7 Conclusions 

The large D expansion of black hole perturbations provides a natural means to isolate the 
sector of quasinormal modes that usually contains the most interesting dynamics of the 

®Of course the hydrodynamic expansion is not convergent itself [15]. 

®In order to put the results of [12] in this form, we use the identity Hm = m~^ -|- Hm-i- 

naive attempt at resumming the highest powers of k in this series gives a result that does not agree 
well with the large-fc behavior in [14]. 
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black hole. This is the decoupled spectrum. In this article we have shown that the large 
D expansion allows efficient calculation of these modes giving unified, analytic expressions 
for all the neutral, static black holes whose solutions are known in closed form,® up to 
relatively high orders in the expansion. 

These results have allowed us to examine the lack of convergence of the series. This 
is not surprising, since the property of decoupling holds to all perturbative orders inl/D 
but is obviously absent at finite values of D. However, it is less clear whether calculations 
for non-decoupling modes are equally affected by these limitations. 
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A Scalar-type perturbations 


For scalar-type perturbations in the ‘Regge-Wheeler’ formulation of (3.1), the minimum 
of the potential is at R ~ D. Although this is within the near-horizon zone, it requires one 
to deal separately with two (overlapping) regions, one with R = 0(1) and another with 
R = 0{D). Although this can be done [3], it is simpler to work instead with a different 
set of gauge invariant variables, X, Y and Z, in terms of which the equations are [13] 


CO 


X'{r) = — — -X(r) + ( Y(r) -h (— - 


LxP' 

2/(r) f[ry 

Z'{r) = X{r) 

together with the consistency condition 

1 


Z{r), 


Z(r), (A.l) 

(A.2) 

(A.3) 


2R 


+ K\Y + — I (D - 2)(D - 3)K - (D - 1)(D - 2)XY - 


2 {D-1){D-3) 


2R 


+ 


- q^f{r) + {D- 2)K'^ - {D - 3)K\y - 


2 2{D-2) , (D-1) 


2n 


R 


+ 


4R2 


Y(r) 


{D - 2)a;V2 “ ^ = 0- 


X{r) 


(A.4) 


“We have not considered static black rings nor other static blackfolds in de Sitter space [16, 17], which 
are not known in closed exact form. 
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(A.5) 


If we introduce two variables P{R) and (5(R) as 

A(R) = P(R) + p:^Q(R), ^(R) = ^R) - 

then the perturbation equations in the near-horizon zone decouple in these variables at 
each order in 1/D. Then they can be solved to yield decoupled quasinormal modes without 
needing to split into two zones. 

B Another perturbative formulation for AdS black branes 

The perturbation equations for AdS black branes can be analyzed in a different manner, 
as described in [14], 

We write the background metric as 

ds2 ^ r^[-f[r)dt^ + dabdx^dx’^) + , (B.l) 

where /(r) = 1 — and the coordinates span For linearized perturbations 

we can always choose a coordinate z aligned with the momentum of the perturbation, 
kaX°' = kz, so the dependence on takes the form ~ Then perturbations can 

be decomposed into scalars, vectors and tensors with respect to the little group of A:“, 
SO{D — 3). In the following we write 

X^ = {z,x^) {i = l,...,D -3). (B.2) 

It is easy to check the equivalence between these perturbations and those of [13]. For 
example, the metric perturbation for a M^“^-scalar perturbation in [14] is, in momentum 
space, 

= ST{t,r)(^kakb-^^^ . (B.3) 

When we align ka with k^ the only non-vanishing component of h^// is hfj, which is 
scalar-type with respect to SO{D — 3). One can find similar relations for other variables 
in the two decompositions. 

The perturbation equations can be written in terms of master variables Zg ■ For tensors 
(s = 2) and vectors (s = 1) these are 

Zs=2'Tij — kij, 

Zg=\di — khii “t" 

while for scalars (s = 0) 

Zg=Q — Aujkhtz “t" 2u^ 

+ (k\{D 3)/(r)) - 2^2) + 2k^htu (B.6) 


(B.4) 

(B.5) 
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Here is a symmetric traceless tensor. Note that since we choose Zg ~ there are no 
scalar-derived tensor or vector perturbations, nor vector-derived tensor perturbations. 
The perturbation equation for tensors Zs =2 is 


7" 

^s=2 


D-1 + f{r) , - k‘^f{r) 

■I-- ^s=2 T-TTTTT;- ^s=2 = U, 


rf{r) 


^f{r) 


(B.7) 


where the prime is derivative with respect to r. 
The equations for the vectors and scalars are 


and 


where 


7" 

^s=\ 


+ 


" + m-f, ] z',_, + ^LzSiZlz.., = 0 


z 


// 

s=0 


+ 


yiA;2 + y2W%, , + Yik^ + 2{D -2)UJ^ 

rf{r)X r^f{r)‘^X 


= 0 , 


X = 2{D- 2)UJ^ - {{D -1) + {D- 3)/(r))A;2, 

= _( 2 Z) _ i)(D - 3)/(r)2 - {D - 1){D - 1 - (D - 4)/(r)), 
Y 2 = 2{D-2){D-l + fir)), 

ys = -{D - 3)r\f\r)ff{r) - {{3D - 7)/(r) + {D - l))u;^ 

Y, = {{D-l) + {D-3)f{r))f{r)- 


(B.8) 


(B.9) 


(B.IO) 


These equations can now be solved perturbatively in 1/D in the near-horizon zone in the 
usual manner, yielding the results in sec. 4.1. 
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C c<;( 3 ) for scalar modes 


Rew(3) 


1 


24q^^Rl 


-Aq^^Rl {2KRo (9 In Rq + 24C(3) + 27r^ + 27) - 3 (SlnRo - 8C(3) + + 2l)) 

+2fRo(^8K^Rl (24C(3) +7r2) 

+12K‘^RI (3 In^ Ro + (21 + 2tt‘^) In Rq + 48C(3) + 2Tr‘^ + 2l) 

-4R:Ro (18 In^ Ro + 3 (47 + 2Tr^) In Rq - 12C(3) + 3l7r2 + 159) 

+48 In^ Ro - 8 { 211 “^ - 45) In Rq - 96C(3) - eOvr^ + 477^ 


-f (^15 + 192i£®R^C(3) + IGK^Ro ((6 + In Ro + 2 (63C(3) + bv^)) 

+2AK^Rl (In^ Ro + 17 In^ Ro + (23 + 27r^) In Ro + 72C(3) - + l) 

-AK'^RI (12 In^ Ro + 195 In^ Ro + 8 (60 + 7^“^) In Ro + 192C(3) + 1847r2 + 57) 


-2R:Ro (-16 In^ Ro - 252 In^ Ro + (l67r2 


969) In Ro + 96C(3) + 527r2 - 



+2K^q^Ro (^16R:^R^ (tt^ In Ro + 30C(3) + 27r^) 

+8R:^R|] (18 In^ Ro + 12 (3 + vr^) In Ro + 198C(3) + 177r2) 

-12il'2R2 (_ in3 Rq + 11 In2 Ro + (63 + IOtt^) In Ro + 8C(3) + 377r2 + 2) 
-4 R:Ro (41n3 Ro + 18 In^ Ro + (227r2 - 135) In Ro + 96C(3) + 6l7r2 - 234) 

+16 In^ Ro + 180 In^ Ro + 501 In Ro + 324^ 


-AK^fRl (^24i£3R3 (in2 Ro + (2 + vr^) In Ro + 2 (8C(3) + vr^)) 

+4^2r 2 (ein^ Ro + (tt^ - 12) InRo - 9 (vr^ - 8C(3))) 

-2KRo (- In^ Ro + 51 In^ Ro + 6 (l3 + bir"^) In Ro + 10 (l2C(3) + - 9)) 

+4 In^ Ro + 99 In^ Ro + 354 In Ro + 234^ 


+8K^Rl i^AK^Rl (3 In^ Ro + 2 (3 + vr^) In Ro + 4 (6C(3) + vr^)) 
—4KRo (6 In^ Ro + 3 (6 + vr^) In Ro + 12C(3) + 87r^) 

+ In^ Ro + 18 In^ Ro + 72 In Ro + 48^ , 


(C.l) 
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Ima;(3) = 4^®RoC(3) 

3 (vr^ In Ro + 8C(3) + - 2KRo (vr^ In Rq + 30C(3) + 27r'^)^ 

+ ^ (^UK^Rl (ln2 Ro + (2 + vr^) In Rq + 2 (8C(3) + vr^)) 

-2KRo (12 In^ Rq + (42 + Ovr^) In Rq + 48C(3) + 297r2) 

+ In^ Ro + 21 In^ Rq + 102 In Ro + 90^ 

-j(^4K^Rl (31n2 Ro + 2 (3 + vr^) In Ro + 4 (6C(3) + vr^)) 

— 4i^Ro (6 In^ Rq + 3 (6 + vr^) In Ro + 12C(3) + 87r^) 

+ ln3Ro + 181n2Ro + 721nRo + 48 ) . (C.2) 



D Results for (A)dS-Schwarzschild black holes 


In this case K = 1 and the eigenvalue must be expanded in 1/n like in (4.7). All the 
dependence on the cosmological constant is included in 


Ro = 


1 

1 - A 


L" 


(D.l) 


where in the last expression we employ A = which is convenient for AdS with 

radius L. 


Vector-type 

wro = —i{i — 1) 

(( 


! + -(£-!-In Ro) 
n 


+ “2 


1 (£_i)_(£_3 + 2Ro)lnRo + ^(lnRo)2 


+ ^ (4(£ - 1) (7r2Ro(2Ro + ^ - 3) - 6Rg(£ - 1)C(3) - 6RoC(3) + 6) 

+2 In Ro (-2 (6 + vr^) R^ + 3 (lO + vr^) Ro + (27r2Ro - 3 (2 vr^) Ro + 12) £ - 24) 
+3(lnRo)2(-4R2 + 8Ro + ^-5)-(lnRo)3)+0(n-^) , (D.2) 
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Scalar-type 


Re uj±ro = ±a/ ^ p 1 + ^ . ((2Ro + 1)£ - 4Ro - 2 Ro In Rq) 

V Rq [ 2 (t - Rojn 

+ 24 (/_Vo) 2 „ 2 (i 6 Ro (^^Ro - 6 ) - 3 (4R2 - 12Ro + l) 

+ (4 (3 + 27r2) r 2 - (132 - Svr^) Rq - 9) f - 4Ro {2tt^RI + 6 (vr^ - 3) Rq - 33) £ 

+ 12Ro In Ro (4(Ro - 2)Ro + (2Ro - 5)£‘^ + (11 - 4Ro)£) - 12 Ro(ln Ro)2(Ro + (Rq - 2)^)) 

+ 48(/_4.()3„3 ((3 {8Ro - 20RS + 6R„ + !)<=- 2«‘{-12 - 15) R„ 

+ 8Rl (12C(3) + 6 + TT^) - 12R2 (-4C(3) + 15 + Stt^) - 3) 

+ £^{-4 (51 + 647r2) Rg + 16R^ (24C(3) + tt^) 

- 24R;] (-40C(3) - 1 + vr^) - 2Ro (96C(3) - 447 + 527r2) -h 15) 

- 8Ro£^(24R^C(3) + 2Rl (l02C(3) + Ttt^) + Rg (96C(3) - 6 - 3l7r2) 

+ Ro (-96C(3) + 114 - eiTT^) 81) 

+ 8 Rg£ ( 8 RD (12C(3) + TT^) + 2R2 (96C(3) + Tvr^) 

-2Ro (60C(3) - 21 + 407r2) + 117) 

- 128R3 ( 7 r 2 (Ro - 2)Ro + 6R2C(3) - 3RoC(3) + 3)) 

- 2Ro In Ro(16R§ (2 (3 + vr^) Rg - 3 (6 + vr^) Rq + 18) + (36Rg - 84Ro + 69) £‘^ 

- 2 ( 47 r 2 (3Rg - 3Ro - 2 ) 3 (24Rg - 66 Ro + 67)) £^ 

+ (8 (6 + bn'^) Rg - 36Rg - 4 (45 + 227r2) Rq + 50l) £‘^ 

- 4Ro ( 47 r 2 R 3 + 2 (24 + Rg - 30 (5 + tt^) Rq + 177) £) 

+ 12 Ro(ln Ro)^(-4Rg ( 2 Rg - 4Ro + 3) + ( 6 Rg - 13Ro + lO) £^ + (-24Rg + 44Ro - 30) £^ 
+ Ro (24R2 - 46Ro + 33) £) - 8 Ro(ln Ro)^ (Rg + (SRg - 6 Ro + 4) £^ - 2Ro£)) 

+0{n-^) , (D.3) 
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and 


Imcjiro = —i{i — 1) 


1 


1 + -(£-2-lnRo) 




1 


+ ^ ' 4 - 3£ + (^ - 2) 


n 

7r2R, 


0 _(2Ro + £-4)lnRo + ^^ 




TT 


TT 


+ ^ { 2Rot — - 2RoC(3) ]+i( 8Ro(2Ro - 1)C(3) + 7 - —Ro(13 - 4Ro) 


TT 


+8 ( y (2 - Ro)Ro - Ro(2Ro - 1)C(3) - 1 


- In Ro ( 4 ( 1 + y ) Rg - 2 (6 + TT^) Ro - 


2Rg7r2 


— (2 + Rq + 5^ £ + 12^ 


+ ^(ln Ro)2 (-4R2 + 8Ro + £ - 6) - + 0{n-^) 


(D.4) 
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